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Abstract 



In this paper we briefly outline the quadrature method for estimating un- 



\0 ! certainties in a function which depends on several variables, and apply it to 

2 ' estimate the numerical uncertainties in QCD-QED rescaling factors. We eni- 

Q ■ ploy here the one-loop order in QED and three-loop order in QCD evolution 
equations of the fermion mass renormalization. Our present calculation is 

Ph| found to be new and also rehable compared to the earher values employed 

JlLi' by various authors. 
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1 Introduction 

An important aspect of particle physics is the concept of mass of a fermion. 
The experimental values of the fermion masses are the input values that we 
always take in various calculations, and they are generally known as physical 
masses of the fermions. Mathematically such physical mass is the running 
mass mf{n) of the fermion defined at the scale equal to its own physical mass, 
i.e., mflu = rrif) = mf{mf). This is true for heavier quarks f = t,b,c which 
have mf{mf) > IGeV and for charged leptons {T,fi,e). However for lighter 
quarks (/ = s, d, u) which have masses lesser than IGeV, the physical mass 
is defined as mf{lGeV). On the other hand, the fermion masses predicted in 
Grand unified theories (GUTs)[l], are defined at the grand unification scale 
Mx- In order to express these masses at lower energy scales we need the 
renormalisation of fermion masses from high scale to low scale[2]. At low en- 
ergy scale below the top-quark mass scale m^, the fermion mass renormalisa- 
tion is governed by the QCD-QED symmetry gauge group SU{?>)c x U{\)em- 
Thus the fermion mass renormalisation from top-quark mass scale down to 
the physical quark mass scale is parametrised by the QCD-QED rescaling 
factors rjf defined by 77/ = mf{mf)/mf{mt). This is an important parame- 
ter which appears in many expressions related to the predictions of GUTs 
at low energy scale. The calcualtion of rjf involves the input values of e.m 
and strong gauge coupling constants a and ^3. The uncertainties associated 
with the experimental values of a and 0:3 will propagate through several in- 
termediate scales upto low energies in the definition of QCD-QED rescaling 
factors, making the uncertainties larger and larger. In this respect, a careful 
estimation of uncertainties in rjf using a reliable method which can control 
the magnifying tendency of the uncertainties in r^j, is highly desirable. The 
uncertainties in rjf so far reported in the literature[4,12], do not match with 
each other, and the methods adopted by them are also not clearly speci- 
fied. In this context we find the quadrature method [3] for the estimation of 
uncertainties of a function which depends on several variables, to be quite 
satisfactory and appropriate. It will be an important numerical exercise to 
estimate the uncertainties in rjf using the quadrature method and compare 
the results with the earlier values. We shall employ the one-loop order in 
QED and three-loop order in QCD evolution equations of the fermion mass 
r enormalizat ion . 

The paper is organised as follows. In Section 2, we outline the procedure 
for the quadrature method. Section 3 is devoted to a brief note on how to 



define QCD-QED rescaling factors and related quantities [4, 5, 6]. In Section 4, 
we present tlie numerical calculations of rescaling factors and their respective 
uncertainties. The paper concludes with a summary in Section 5. 

2 Estimation of uncertainties in Quadrature 
Method 

We define the quadrature method in the following way. If a function F 
depends on several variables Xi, where i=l,2,...n, then the uncertainties in 
F resulting from the uncertainties of the independent variables Xi can be 
estimated by the following expression, 



6F = ± 



\ 



^O''^-'^ 



(1) 



where ± ( ^j {5x,i) is the error in F due to the error in i^^ variable Xi. Here 
the uncertainties (or random errors of independent variables) add in quadra- 
ture[3] like orthogonal vectors, a and b whose resultant magnitude is given 
by vaM-&^. 

We give here three simple properties which will be needed to estimate the 
uncertainties in different quantities leading to the calculations of rescaling 
factors and also other parameters which depend on it (e.g., neutrino masses 
and mixings): 

(i): If F{x,y) = f{x,y) x g{x,y), the uncertainties in F due to the random 
errors ±6x and ±6y in the independent variables x and y respectively, are 
given by 



6F = ± 



\ 
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(ii): If F(x,y) = (^), the uncertainties in F due to the random errors 
±6x and ±6y in the independent variables x and y respectively, are given by 



6F = ± 



\ 
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(l) , then the uncertainties in F due to the un- 
certainties ±6x and ±6y in the independent variables x and y respectively, 



(iii): li F{x,y) 



are 
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3 QCD-QED rescaling factors 



The QCD-QED rescaling factor rjf of the fermion f , which can take care of 
the fermion-mass-renormalisation from the top-quark mass scale down to the 
physical fermion mass scale mf{mf), is defined as [7] 
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for f = b, c quarks having mf{mf) > 1 GeV, 
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for f = s, d, u quarks having mf{mf) < 1 GeV, and 
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Here mf{mt) is the running mass of the fermion f at the top-quark mass scale, 
fi = rrit, and mf{mf) corresponds to the physical mass of the fermion. For 
convenience we also define a quantity R^{fi,fi') which represents the inverse 
of the QCD-QED rescaling factor in the narrow range, /x' — /i, (/i > fj,') , 
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= R-QEDif^, f^') X RQCDif^, f^') (8) 

The definitions of r]f and R-^ in Eqs. (5)-(8), require only the contributions 
from the QED part for the charged leptons. Making use of Eq.(8) for suc- 
cessive narrow mass ranges: mb-nit, nir-mb, rric-mr, IGeV-mc, ms-lGeV, 



m^-nis, 771,1-171^, 171^-171^, 771^-171^, the QCD-QED rescaling factor rjf defined 
in Eqs.(5)-(7), can be rewritten as 

Vb = r]b{mb,mt) = 1 / R'' {mb, rut) , 

Tjr = r]^{m^, rut) = 1/ [RqEDi'^r, mb)RQEDi'mb, mt)\ ,1 = T,fi,e, 
Vc = Vc{mc,mt) = 1/ [R^imc, m^)R''{m^, mb)R''{mb, mt)] , 
Vs,d = r]s,d{lGeV,mt) = r]^/ [R'''^{lGeV,m^)R'^'^{m^,m^)R'^'^{m^,mb) 
Vu = VuilGeV, mt) = r^J [R^{lGeV, mj] , 

^M = Vf^{mfi,m) = Vr/ [RQEDi^f^^^s)RQEDi^s, IGcV)] , 

r]e = r]e{me, mt) = r]^,/ RQE^{me, mu)RQEDi^u, mdjRQEoi^d, rrif,) (9) 
We use one-loop order in QED and three-loop order in QCD evolution equa- 
tions of fermion mass renormalistion for the evaluations of Rqed ^^^ ^qcd 
respectively. The contribution of one-loop QED, running from the scale n' 
to fi, to the rescaling factors through Eq.(8) is now given by [7] 



where 
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(11) 

The summation in Eqs.(ll) is over the active fermions at the relevant mass 
scale, and f is the specific fermion under consideration. We employ the one 
loop RGE for the estimation of e.m. gauge couplings Q;(/i) at successive 
renormalisation points. 
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The three-loop QCD running quark mass formula is given by [7,^ 
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where m„ is the scale-invariant quark mass, and 
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The QCD /3 -functions and anomalous dimensions are given as, 

70 = 4, 
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where n/ is the number of quark flavours at the relevant mass scale and 
^(3) = 1.202. The QCD rescaling contribution to R^ in the relevant mass- 
range (fi'-fJ,) defined in Eq.(8), can be obtained from Eq.(13) as 
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since 0:3 (/x) changes smoothly with //, so values of the constants A and B 
remain as same for a given energy range including the limits. 

Using Eqs.(lO) and (14) via (8), the QCD-QED rescaling factors can be 
estimated. The values of a^{^) in Eq.(16) can be obtained by solving the 
three-loop QCD RGE with (ag = ^3/47r),[9] 
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The conventional solution of Eq.(17) having a constant of integration called 
QCD dimensional parameter A-which provides a parametrisation of the n 
dependence of 0:3 (/i) [9] is given by 
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where 



Since /3-function coefficients change by discrete amount as flavour thresholds 
are crossed while integrating the differential equation(17) for a^i^fi), A also 
changes to take care of the validity of Eq.(18) for all values of fi. This leads 
to the concept of different A*^"^-'') for each range of n corresponding to an 
effective number of quark flavours Uf. In the MS scheme, one finds the 
relation among different A*^"-^) [10] as 
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The QCD parameter A)Jg which corresponds to n/ = 5 can be calculated by 
using Eq.(18) with the input experimental value of as{Mz) and then A^^ 



,(3) 



and Ajjg can be obtained from Eq.(20) through a search programme. Sub- 
sequently, 0:3 (/i) and Rqcd can be calculated for the energy, n < Mz- 



If we consider SUSY above Mz-, assuming the existence of one-light Higgs 
doublet {Nh) and five quark flavours (nj = 5) in the energy range irtt — Mz-, 
the strong gauge coupling at the scale rrit are evaluated using RGE solution 
with one loop order 
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for li > Mz- 



4 Numerical calculations and results: 

For numerical calculation we take the input values of fermion masses as [9] 

mt = l75GeV,mb = 4.25^61/, m^ = l.785GeV, 

rric = 1.25Ge1/,m, = 0.18GeV,m^ = 0.l05GeV, 

rrid = 0.009GeV, m„ = O.OObGeV, m^ = .OOOblGeV, (22) 

and CERN-LEP data [11] of a-^{Mz) = 127.9 ± 0.1 and strong coupling 
constant aaiMz) = 0.1172 ± .002 at Mz = 91.187G'e\/, which is referred to 
as Case-I. For comparison, we also estimate the uncertainties using values 
of a3{Mz) = 0.120 ± .0028 referred to as Case-II and a3{Mz) = .118 ± .007 
made use for earlier calculation of uncertainties [12], now referred to as Case- 
III. 

Now using above data, we evaluate the gauge couplings a{fi) and 0:3 (/x) 
at various renormalisation points starting from rrit down to the individ- 
ual fermion mass (for quark it stops at IGeV) from their corresponding 
Eqs.(12),(18)and (21) and values are presented in Table-1 and Table-5 re- 
spectively. The co-efficients of /9-functions and anomalous dimensions in 
RGEs for QED and QCD with constants A, B relevant in different energy 
ranges are estimated and shown in Tables- 1,2, 3. We also calculate the in- 
verse of the fermion mass renormalisation factors Rf {jj,, fi') for charged 

leptons, up-quarks and down-quarks shown in Table-1, 2 and Ri {l^^ 1^') in 
three Cases-I, II, III in Table-5. The values of QCD dimensional parameters 
^^Ms fo^ estimating 0:3 (/x) at various energy scales are presented in Table-4. 
We present below the numerical values of QCD-QED rescaling factors rjf 
(f=b, c, s, d,u) with uncertainties for different values of a^i^Mz). It is ob- 
served that the uncertainties in the input values of 0:3 (M^) propagate down 
to the lower energy scale, causing the uncertainties in rjj larger and larger: 



Case asjMz) rjb Vc Vs^ Vu 

I 0.1172 ±.002 1.5380t[J:|]ii 2.1785t''oilm 2.4329l[iJ?^i 2.4489lg|?^^ 



II 0.1200 ±.0028 1.5649lH^?^ 2.2487l[l:ii? 2.6486t[l:^iJ 2.6660t[J:il 

III 0.1180 ±.007 1.5458l°:l?i 2.2110l°:118 2.5106lJ:^?^^ 2.527llJ:^?^^ 



The rescaling factors for charged leptons are found as 

7]r = 1.017 ± .0007, T]^ = 1.027 ± .0038, r]e = 1.046 ± .0099 for all cases. 

For Case-Ill with the value of q;3(M^) = .118 ± .007, the results estimated 
by Deshpande and Keith [12] using the same value of 0:3 (M^) = .118 ± .007, 
have small uncertainties as follows. 

f], = 1.53iH^, Vc = 2.09t°:S, Vs,ci = 2.36l[J:i, Vu = 2.38l[JJ^. 

These low uncertainties may be resulted from an apparent mistake in 
their estimation of uncertainties in 
o^s^if^t) = 9.30+0:05!, 
instead, it should be corrected as 

a^\fit) = 9.30;°:tl. 

which will lead to larger uncertainties in r]f. 

The present estimation of uncertainties in r]f using quadrature method is 
very convincing and also it regulates the propagation of uncertainties while 
running the energy scale from high to low scales. These uncertainties can 
reliably be used in other low energy calculations e.g., neutrino masses and 
mixings in see-saw mechanism in low enery scale [13]. 



5 Summary and conclusion: 

To summarise, we have outlined the procedure for estimating the uncertain- 
ties using the quadrature method. In particular, we employ this method to 



estimate the numerical uncertainties in QCD-QED rescaling factors rjf ac- 
quired from the uncertainties in input values of a^lMz) and a{Mz) while 
running the energy scale from high to low. We have used the three-loop or- 
der in QCD and one-loop order in QED evolution equations to calculate the 
rescaling factors while their uncertainties estimated by the qudrature method 
are found to be new in comparison to earlier estimates. 
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Table-1: Coefficients of /9-functions in the RGEs for QED,the values of 
Tj in different energy ranges (/i — /i'), /i > /i', the inverse of the gauge 
couphng a~^(yu) and inverse of the rescahng factors Ri {l^^l^') for charged 
leptons (e,/i,r) with input value of a~^{Mz) = 127.9 ± 0.1. 



Energy range 




(/ = e,/i,r) 


a-^/i) 


^QED 

(/ = e,/i,r) 


mt — rUb 


80/9 


-0.3375 


126.98 ±0.1 


0.9865 ± 0.0002 


iTib — rrir 


79/9 


-0.3553 


132.24 ±0.1 


0.9968 ± 0.0005 


rrir — rrir 


64/9 


-0.4219 


133.40 ±0.1 


0.9987 ±0.0006 



mr-lGeV 16/3 -0.5625 133.81 ±0.1 0.9990 ± 0.0008 

IGeV-ms 16/3 -0.5625 133.99 ±0.1 0.9939 ± 0.0007 

rris-m^ 44/9 -0.6136 135.45 ±0.1 0.9981 ± 0.0089 

m^-rrid 32/9 -0.8438 135.87 ±0.1 0.9914 ± 0.0011 

md-m^ 29/9 -0.9643 137.26 ±0.1 0.9979 ± 0.0014 

rriu-me 4/3 -2.2500 137.55 ±0.1 0.9921 ± 0.0032 

m, -0 ... ... 138.04 ±0.1 
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Table-2: Coefficients of /9- functions in the RGEs for QED,the values of 
Tj in different energy ranges {fi — ji'), /i > /i', and inverse of the rescal- 
ing factors ni {fi, ji') for up- and down-quarks with renormahsation scale 
stopping at IGeV. 



Energy scale 6^" ^ ^ B^d B^d 

(/i — /i') (f=u,c.t) (f=s,d,b) (f=u,ct) (f=s,d,b) 

mt-mi, 80/9 -0.15 -0.0375 0.9839 ± 0.0002 0.9985 ± 0.0001 

mi,-mr 79/9 -0.1579 -0.0395 0.9986 ± 0.0002 0.9997 ± 0.0001 

m^-ruc 64/9 -0.1875 -0.0469 0.9994 ± 0.0003 0.9999 ± 0.0001 

mc-lGeV 16/3 -0.25 -0.0625 0.9996 ± 0.0004 0.9999 ± 0.0001 



Table-3: Coefficients of f3 functions in RGEs, anomalous dimensions for 
QCD and values of the constants A and B defined in Eq.(lO) with 70 = 4 
throughout the computation. 



Energy scale 71 72 &o &i &2 A B nf 
(/^ - ^^') 

mt-mi, 506/9 474.89 23/3 116/3 9769/54 4.7020 24.0123 5 

mi,-mr 526/9 636.62 25/3 154/3 21943/54 4.0563 22.2280 4 

m^ - rric 526/9 636.62 25/3 154/3 21943/54 4.0563 22.2280 4 

mc - IGeV 546/9 794.90 9 64 34767/54 3.5802 21.9434 3 
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Table-4: The values of QCD dimensional parameters A^^^ with uncer- 
tainties estimated for three different values of 0:3 (/u) for cases I, II, III. 



case AMgGeV AmsGcV AmsGcV 

I 0.1995tnif, 0.2877lHii 0.3375i°:°ii 

TT n 9Q9Q+0.0364 r^ qqnc:+0.0459 r, qoqo+0.0472 

ii U./:d/y_o.o334 U.ddUO_o.o429 U.d8d5_o.o448 

III 0.20871°:°??^ 0.2995tr,lll 0M98lliZ 

Table-5: The values of Q;3(/i) and Rq(jj^{^ — ^') with uncertainties es- 
timated for quarks using three-loop order in QCD evolution equations for 
three different input experimental values of a^{Mz)- Here QM stands for 
the quadrature method employed here. 

case /parameter mt — rrih rrih — 171^- nij. — rric rric — IGeV IGeV — 

I/a3(/i) 0.10681°:^ 0.2225^^^?^ 0.3149^°;°}^° 0.39671°;°^^ 0.4819l°S 

jQM)/ Rjco 0.6512±g:°li 0.81851°:°^^° 0.8682±°:°ig 0.8900±°:g^°° 

II/a3(/i) 0.109liHgi 0.2334iHl^? 0.33821^11 0.42331°;°^^ 0.55681°;^?°^ 

{QM)/R(,an 0.64001^?^^ 0.80551^1^ 0.86961°:°!^ 0.8439l°:S^l 

m/asif.) 0.l075l:',fil 0.22571°:°!^ 0.32131°:°^^! 0A079t',fit 0.500lt'oitll 
{QM)/Ria^ OM79trolll 0-S152tUlll 0-8633l°:^°I° 0.8753l°:f,i 
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